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The thermodynamic entropy of an isolated system is given by its von Neumann entropy. Over the 
last few years, there is an intense activity to understand thermodynamic entropy from the principles 
of quantum mechanics. More specifically, is there a relation between the (von Neumann) entropy of 
entanglement between a system and some (separate) environment is related to the thermodynamic 
entropy? It is difficult to obtain the relation for many body systems, hence, most of the work 
in the literature has focused on small number systems. In this work, we consider black-holes — 
that are simple yet macroscopic systems — and show that a direct connection could not be made 
between the entropy of entanglement and the Hawking temperature. In this work, within the 
adiabatic approximation, we explicitly show that the Hawking temperature is indeed given by the 
rate of change of the entropy of entanglement across a black hole’s horizon with regard to the 
system energy. This is yet another numerical evidence to understand the key features of black hole 


thermodynamics from the viewpoint of quantum 
PACS numbers: 03.67.Mn. ll.10.-z, 05.50.+q, 05.70.- 

I. INTRODUCTION 

Equilibrium statistical mechanics allows a successful 
description of the thermodynamic properties of matter 
W- More importantly, it relates entropy, a phenomeno¬ 
logical quantity in thermodynamics, to the volume of a 
certain region in phase space [5]. The laws of thermody¬ 
namics are also equally applicable to quantum mechan¬ 
ical systems. A lot of progress has been made recently 
in studying the cold trap atoms that are largely isolated 
from surroundings EH5]. Furthermore, the availability 
of Feshbacli resonances is shown to be useful to control 
the strength of interactions, to realize strongly correlated 
systems, and to drive these systems between different 
quantum phases in controlled manner |TT)TfTm . These 
experiments have raised the possibility of understand¬ 
ing the emergence of thermodynamics from principles of 
quantum mechanics. The fundamental questions that 
one hopes to answer from these investigations are: How 
the macroscopic laws of thermodynamics emerge from 
the reversible quantum dynamics? How to understand 
the thermalization of a closed quantum systems? What 
are the relations between information, thermodynamics 
and quantum mechanics [I4l - [T9] ? While answer to these 
questions, for many body system is out of sight, some 
important progress has been made by considering simple 
lattice systems (See, for instance, Refs. [20l - (23] ). In this 
work, in an attempt to address some of the above ques¬ 
tions, our focus is on another simple, yet, macroscopic 
system — black-holes. 

It has long been conjectured that a black hole’s ther¬ 
modynamic entropy is given by its entropy of entangle¬ 
ment across the horizon [MKSIj- However, this has never 
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information theory. 


been directly related to the Hawking temperature mi¬ 
ller e we show that: 

(i) Hawking temperature is given by the rate of change 
of the entropy of entanglement across a black hole’s 
horizon with regard to the system energy. 

(ii) The information lost across the horizon is related to 
black hole entropy and laws of black hole mechanics 
emerge from entanglement across the horizon. 

The model we consider is complementary to other 
models that investigate the emergence of thermodynam¬ 
ics HHIS]: First, we evaluate the entanglement entropy 
for a relativistic free scalar fields propagating in the 
black-hole background while the simple lattice models 
that were considered are non-relativistic. Second, quan¬ 
tum entanglement can be unambiguously quantified only 
for bipartite systems [32J [33]. While the bipartite system 
is an approximation for applications to many body sys¬ 
tems, here, the event horizon provides a natural bound¬ 
ary. 

Evaluation of the entanglement of a relativistic free 
scalar field, as always, is the simplest model. However, 
even for free fields it is difficult to obtain the entangle¬ 
ment entropy. The free fields are Gaussian and these 
states are entirely characterized by the covariance matrix. 
It is generally difficult to handle covariance matrices in 
an infinite dimensional Hilbert space [33 :. There are two 
ways to calculate entanglement entropy in the literature. 
One approach is to use the replica trick which rests on 
evaluating the partition function on an n-fold cover of the 
background geometry where a cut is introduced through¬ 
out the exterior of the entangling surface |33l |34]- Second 
is a direct approach , where the Hamiltonian of the field is 
discretized and the reduced density matrix is evaluated in 
the real space. We adopt this approach as entanglement 
entropy may have more symmetries than the Lagrangian 
of the system [351 . 
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To remove the spurious effects due to the coordinate 
singularity at the horizon 1 , we consider Lemaitre coordi¬ 
nate which is explicitly time-dependent [57] . One of the 
features that we exploit in our computation is that for a 
fixed Lemaitre time coordinate, Hamiltonian of the scalar 
field in Schwarzschild space-time reduces to the scalar 
field Hamiltonian in flat space-time [25] . 

The procedure we adopt is the following: 

(i) We perturbatively evolve the Hamiltonian about 
the fixed Lemaitre time. 


(ii) We obtain the entanglement entropy at different 
times. We show that at all times, the entanglement 
entropy satisfies the area law i. e. 5(e) = C(e)A 
where 5(e) is the entanglement entropy evaluated 
at a given Lemaitre time (e), C(e) is the propor¬ 
tionality constant that depends on e, and A is the 
area of black hole horizon. In other words, the value 
of the entropy is different at different times. 

(iii) We calculate the change in entropy as function of 
e, i. e., AS/Ae. Similarly we calculate change in 
energy E(e), i.e., AE/Ae. 


For several black-hole metrics, we explicitly show that 
ratio of the rate of change of energy and the rate of change 
of entropy is identical to the Hawking temperature. 

The outline of the paper is as follows: In Sec. @, 
we set up our model Hamiltonian to obtain the entangle¬ 
ment entropy in (D + 2)-dimensional space time. Also, 
we define entanglement temperature , which had the same 
structure from the statistical mechanics, that is, ratio of 
change in total energy to change in entanglement entropy. 
In Sec. (III), we numerically show that for different black 
hole space times, the divergent free entanglement tem¬ 
perature matches approximately with the Hawking tem¬ 
perature obtained from general theory of relativity and 
its Lovelock generalization. This provides a strong ev¬ 
idence towards the interpretation of entanglement en¬ 
tropy as the Bekenstein-Hawking entropy. Finally in Sec. 
(IVI, we conclude with a discussion to connect our analy¬ 
sis with the eigenstate thermalization hypothesis for the 
closed quantum systems [221 . 

Throughout this work, the metric signature we adopt 
is (+, —, —, —) and set h = ks = c = 1. 


II. MODEL AND SETUP 
A. Motivation 

Before we go on to evaluating entanglement entropy 
(EE) of a quantum scalar field propagating in black-hole 
background, we briefly discuss the motivation for the 


1 In Schwarzschild coordinate, r > 2 M need to be bipartited m- 


studying entanglement entropy of a scalar field. Consider 
the Einstein-Hilbert action with a positive cosmological 
constant (|A|): 


S EH (g) = Ml I d A x^=g [R - 2|A|] . (1) 


where R is the Ricci scalar and M P1 is the Plank mass. 
Perturbing the above action w.r.t. the metric g= g^-\- 
h^, the action up to second order becomes [28]: 

M 2 I" 

S EH (g,h) = —f-J d A xs/\g\ [V a VV a /i F 

+ |A|V^H- (2) 

The above action corresponds to massive (A) spin-2 field 
(h^v) propagating in the background metric g Rewrit¬ 
ing, [where e M „ is the constant po¬ 

larization tensor], the above action can be written as 

S EH (9,h) =~\J d 4 xV\g\ + |A|4> 2 ] . (3) 

which is the action for the massive scalar field propa¬ 
gating in the background metric g In this work, we 
consider massless (A = 0 corresponding to asymptoti¬ 
cally flat space-time) scalar field propagating in (D + 
2)—dimensional spherically symmetric space-time. 


B. Model 

The canonical action for the massless, real scalar field 
$(a: Al ) propagating in (D + 2)—dimensional space-time is 

s=\j d D+2 x yf^g^d^x) 3 v $(x) (4) 

where g^ is the spherically symmetric Lemaitre line- 
element m- 

ds 2 = dr 2 - (1 - /[r(r,£)]) d£ 2 - r 2 {r,t)dtt 2 D (5) 

where r, £ are the time and radial components in 
Lemaitre coordinates, respectively, r is the radial dis¬ 
tance in Schwarzschild coordinate and dflo is the 
D —dimensional angular line-element. In order for the 
line-element <[5j) to describe a black hole, the space-time 
must contain a singularity (say at r = 0) and have hori¬ 
zons. We assume that the asymptotically flat space-time 
contains one non-degenerate event-horizon at r^. The 
specific form of f(r) corresponds to different space-time. 

Lemaitre coordinate system has the following interest¬ 
ing properties: 

(i) The coordinate r is time-like all across 0 < r < oo, 
similarly £ is space-like all across 0 < r < oo. 

(ii) Lemaitre coordinate system does not have coordi¬ 
nate singularity at the horizon. 
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(iii) This coordinate system is time-dependent. The test 
particles at rest relative to the reference system are 
particles moving freely in the given field m- 

(iv) Scalar field propagating in this coordinate system is 
explicitly time-dependent. 

The spherical symmetry of the line-element ([5]) allows 
us to decompose the normal modes of the scalar field as: 


C. Important observations 


$ ( X ) = ^2 Z lmi {0,<t>i), 


( 6 ) 


where i £ {1, 2,... D — 1} and Z lm ’s are the real hyper- 
spherical harmonics. We define the following dimension¬ 
less parameters: f = r/r h , I = t;/ r h, f = r/r h , &i m = 
rh^im- By the substitution of the orthogonal proper¬ 
ties of Z lm _, the canonical massless scalar field action 
becomes, 


5 = f dr d£\ 


: D 




1 


J 2 - V 1 ~ W] l(yl+ p - ~'^Irr 


v 71 - /[ f ] 

(7) 


The above action contains non-linear time-dependent 
terms through f(r). Hence, the Hamiltonian ob¬ 
tained from the above action will have non-linear time- 
dependence. While the full non-linear time-dependence 
is necessary to understand the small size black-holes, 
for large size black-holes, it is sufficient to linearize the 
above action by fixing the time-slice and performing the 
following infinitesimal transformation about a particular 
Lemaitre time f [5H]. More specifically, 


t —t t = r + e, 

-/ r f\ 


£ — ► — C, (8a) 

= r(f + e,£), (8b) 

(r, £) -> $\ m . {?, £') = §i mi (?, £) (8c) 

where e is the infinitesimal Lemaitre time. The functional 
expansion of /(f) about e and the following relation be¬ 
tween the Lemaitre coordinates m, 


€-f= 


dr 


\J 1 - f[r{f,£)] 


(9) 


The Hamiltonian in Eq. (10) is key equation regarding 


which we would like to stress the following points: First, 
in the limit of e —> 0, the Hamiltonian reduces to that 
of a free scalar held propagating in hat space-time [28] , 
In other words, the zeroth order Hamiltonian is identical 
for all the space-times. Higher order e terms contain 
information about the global space-time structure and, 
more importantly, the horizon properties. 

Second, the Lemaitre coordinate is intrinsically time- 
dependent; the e expansion of the Hamiltonian corre¬ 
sponds to the perturbation about the Lemaitre time. 
Here, we assume that the Hamiltonian H undergoes adi¬ 
abatic evolution and the ground state is the instan¬ 
taneous ground state at all Lemaitre times. This assump¬ 
tion is valid for large black-holes as Hawking evaporation 
is not significant. Also, since the line-element is time- 
asymmetric, the vacuum state is Unrulr vacuum. Eval¬ 
uation of the entanglement entropy for different values 
of e corresponds to different values of Lemaitre time. As 
we will show explicitly in the next section, entanglement 
entropy at a given e satishes the area law [5(e) oc A] 
and the proportionality constant depends on e i. e. 
5(e) = C(e)A. 

Third, it is not possible to obtain a closed form analytic 
expression for the density matrix (tracing out the quan¬ 
tum degrees of freedom associated with the scalar field in¬ 
side a spherical region of radius 77 ,) and hence, we need to 
resort to numerical methods. In order to do that we take 
a spatially uniform radial grid, {rj}, with b = rj + 1 — Vj. 
We discretize the Hamiltonian H in Eq. (101. The proce¬ 
dure to obtain the entanglement entropy for different e is 
similar to the one discussed in Refs. ;25] [28]. In this 
work, we assume that the quantum state correspond¬ 
ing to the discretized Hamiltonian is the ground state 
with wave-function fl/GS^i!..., x n ; j/i,..., j/jv-n)- The 
reduced density matrix p(y,y') is obtained by tracing 
over the first n of the N oscillators 

/ n 

(J| dxt) ^gs(x 1 , x n ;y)^* GS {xi ,.., x n \y) 

i—1 

(u) 

Fourth, in this work, we use von Neumann entropy 

5(e) = —TV (p log p) (12) 


allow us to perform the perturbative expansion in the 
above action. 

After doing the Legendre transformation, the Hamil¬ 
tonian up to second order in e is 

H(e) ~H 0 + eV 1 + e 2 V 2 (10) 

where H 0 is the unperturbed scalar field Hamiltonian in 
the flat space-time, V 1 and V 2 are the perturbed parts of 
the Hamiltonian (for details, see Appendix |A|. Physi¬ 
cally, the above infinitesimal transformations T8q) corre¬ 
spond to perturbatively expanding the scalar field about 
a particular Lemaitre time. 


as the measure of entanglement. In analogy with mi- 
crocanonical ensemble picture of equilibrium statistical 
mechanics, evaluation of the Hamiltonian H at different 
infinitesimal Lemaitre time e, corresponds to setting the 
system at different internal energies. In analogy we define 
entanglement temperature [39] : 

1 _ A5(e) Slope of EE(A5/Ae) 

T ee A E{e) Slope of energy (AE/ Ae) 

The above definition is consistent with the statistical 
mechanical definition of temperature. In statistical me¬ 
chanics, temperature is obtained by evaluating change in 
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the entropy and energy w.r.t. thermodynamic quantities. 
In our case, entanglement entropy and energy depend 
on the Lemaitre time, we have evaluated the change in 
the entanglement entropy and energy w.r.t. e. In other 
words, we calculate the change in the ground state en¬ 
ergy (entanglement entropy) for different values of e and 
find the ratio of the change in the ground state energy 
and change in the EE. As we will show in the next sec¬ 
tion, EE and energy goes linearly with e and hence, the 
temperature does not depend on e. While the EE and 
the energy diverge, their ratio is a non-divergent quan¬ 
tity. To understand this, let us do a dimensional analysis 


/ T \ D+l A 

[Ed+i] oc N d+1 oc , [S'] a (14a) 

=> [Tee] oc ^ oc N — (14b) 

= [Tee] oc ( N/n) D => finite (14c) 


where Ap is the D+l dimensional hyper- surface area. In 
the thermodynamic limit, by setting L finite with N —> 
oo and b — > 0, Tee in Eq. (14c) is finite and independent 
of e. 

For large N, we show that, in the natural units, the 
above calculated temperature is identical to Hawking 
temperature for the corresponding black-hole m-- 


= 1 df(r) 
2n 4-7T dr 


(15) 


Fifth, it is important to note the above entanglement 
temperature is non-zero only for f(r) ^ 1. In the case of 
flat space-time, our analysis shows that the entanglement 
temperature vanishes, and we obtain Tee numerically for 
different black hole space-times. 


III. RESULTS AND DISCUSSIONS 




Figure 1: The plots of total energy (left) and EE (right) as a 
function of e for the 4d Schwarzschild black hole. We set N = 
600 and n = 150. The cyan coloured dots are the numerical 
data and the red line is the best linear fit to the data. 


A. Schwarzschild (SBH) black holes 


The 4-dimensional Schwarzschild black hole space-time 
( put D = 2) in dimensionless units f is given by the line 
element in Eq. !> with /(f) is given by: 


/(f) = 1 - \ 


(16) 


In Fig.Q, we have plotted total energy (in di¬ 
mensionless units) and EE versus e for 4-dimensional 
Schwarzschild space-time. Following points are impor¬ 
tant to note regarding the numerical results: First for 
every e, von Neumann entropy scales approximately as 
S ~ ( Th/b ) 2 . Second, EE and the total energy increases 
with e. 


Using relation (131, we evaluate “entanglement” tem¬ 


perature numerically. In dimensionless units, we get 
Tee = 0.0793 which is close to the value of the Hawk¬ 
ing temperature 0.079. However, it is important to note 
that for different values of N, we obtain approximately 
the same value of entropy. The results are tabulated, see 
TableQ. See Appendix |C| for plots of energy and EE 
for n = 50,80,100 and 130. 


The Hamiltonian H in Eq. (10) is mapped to a system 
of N coupled time independent harmonic oscillators (HO) 
with non-periodic boundary conditions. The interaction 
matrix elements of the Hamiltonian can be found in Ref 
001. The total internal energy (E) and the entanglement 
entropy (S') for the ground state of the HO’s is computed 
numerically as a function of e by using central difference 
scheme (see Appendix |B|). All the computations are done 
using MATLAB R2012a for the lattice size N = 600, 
10 < n < 500 with a minimum accuracy of 10 -8 and a 
maximum accuracy of 1CU 12 . 

In the following subsections, we compute Tee numer¬ 
ically for two different black-hole space-times, namely, 
4 dimensional Schwarzschild and Reissner-Nordstrom 
black holes and show that they match with Hawking tem¬ 
perature Teh- Tee is calculated by taking the average 
of entanglement temperature for each n’s by fixing N. 


B. Reissner-Nordstrom (RN) black holes 


The 4-dimensional Reissner-Nordstrdm black hole is 
given by the line element in Eq. <©■ where /(f) is 

m = i- 2 +ln + (Qpl (17 , 

Q is the charge of the black hole. Note that we have 
rescaled the radius w.r.t the outer horizon ( 77 , = M + 
yjM 2 - Q 2 ). Choosing q = Q/r h , we get 


f{r) = 1 - 


U + <7 2 ) 


(18) 


and the black hole temperature in the unit of r^ is Tbh = 
(1 - q 2 )/ 4tt. 
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Figure 2: The plots of total energy (top panel) and EE (bottom panel) for as a function of e for different q’s in 4- dimensional 
Reissner-Nordstrom black hole. For these, we set N = 600 and n = 150. The cyan coloured dots are the numerical data and 
the red line is the best linear fit to the data. 


Black hole 
space time 


Tbh 

Tee 

4d-SBH 


0.07958 

0.07927 


9 

= 0.1 

0.07878 

0.07836 

4d-RN 

9 

= 0.2 

0.07639 

0.07507 


9 

= 0.3 

0.07242 

0.07501 


9 

= 0.4 

0.06685 

0.06659 


Table 1: The table lists entanglement temperature and 

Hawking temperature (measure in units of rn ) for 4- 
dimensional Schwarzschild and Reissner-Nordstrom black 
hole space times. 


Note that we have evaluated the entanglement temper¬ 
ature by fixing the charge q. For a fixed charge q , the first 
law of black hole mechanics is given by dE = (k/2tt) dA , 
where A is the area of the black hole horizon. The en¬ 
ergy and EE for different q values have the same profile, 
which looks exactly like in the previous case and is shown 
in the middle row in Fig. ©• See Appendix |C| for plots 
for other values of n. As given in the table (1), Tee 
matches with Hawking temperature. 


IV. CONCLUSIONS AND OUTLOOK 

In this work, we have given another proof that 4- 
dimensional black hole entropy can be associated to en¬ 
tropy of entanglement across the horizon by explicitly 
deriving entanglement temperature. Entanglement tem¬ 


perature is given by the rate of change of the entropy of 
entanglement across a black hole’s horizon with regard 
to the system energy. Our new result sheds the light on 
the interpretation of temperature from entanglement as 
the Hawking temperature, one more step to understand 
the black hole thermodynamics from the quantum infor¬ 
mation theory platform. 

Some of the key features of our analysis are: First, 
while entanglement and energy diverge in the limit of 
6 —» 0, the entanglement temperature is a finite quan¬ 
tity. Second, entanglement temperature vanishes for the 
flat space-time. While the evaluation of the entangle¬ 
ment entropy does not distinguish between the black-hole 
space-time and flat space-time, entanglement tempera¬ 
ture distinguishes the two space-times. 

Our analysis also shows that the entanglement entropy 
satisfies all the properties of the black-hole entropy. First, 
like the black hole entropy, the entanglement entropy in¬ 
creases and never decreases. Second, the entanglement 
entropy and the temperature satisfies the first law of 
black-hole mechanics dE = Tee dS. We have shown this 
explicitly for Schwarzschild black-hole and for Reissner- 
Norstrom black-hole . 

It is quite remarkable that in higher dimensional space 
time the Renyi entropy provides a convergent alternative 
to the measure of entanglement [30], however, entangle¬ 
ment temperature will depend on the Renyi parameter. 
While a physical understanding of the Renyi parameter 
has emerged mi, it is still not clear how to fix the Renyi 
parameter from first principles [42] . 
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Our analysis throws some light on the emergent gravity 
paradigm |43ft46l where gravity is viewed not as a funda¬ 
mental force. Here we have shown that the information 
lost across the horizon is related to the black-hole en¬ 
tropy and the laws of black-hole mechanics emerge from 
the entanglement across the horizon. Since General Rel¬ 
ativity reduces gravity to an effect of the curvature of 
the space-time, it is thought that the microscopic con¬ 
stituents would be the atoms of the space-time itself. Our 
analysis shows that entanglement across horizons can be 
used as building blocks of space-time mm- 

One of the unsettling questions in theoretical physics 
is whether due to Hawking temperature the black-hole 
has performed a non-unitary transformation on the state 
of the system aka information loss problem. Our analy¬ 
sis here does not address this for two reasons: (i) Here, 
we have fixed the radius of the horizon at all times and 
evaluated the change in the entropy while to address the 
information loss we need to look at changing horizon ra¬ 
dius. (ii) Here, we have used perturbative Hamiltonian, 
and hence, this analysis fails as the black-hole size shrinks 
to half-its-size [|5]. We hope to report this in future. 

While the unitary quantum time-evolution is reversible 


and retains all information about the initial state, we 
have shown that the restriction of the degrees of freedom 
outside the event-horizon at all times leads to temper¬ 
ature analogous to Hawking temperature. Our analysis 
may have relevance to the eigenstate thermalization hy¬ 
pothesis [201423] . which we plan to explore. 
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Appendix A: Calculation of Scalar field Hamiltonian in Lemaitre coordinate 

In this Appendix section, we give details of the derivation of the Hamiltonian (H) upto second order in e. Using 
the orthogonal properties of the real spherical harmonics Z lm , the scalar field action reduces to, 


1 


s = ~Y, dfd Z 


~D -rru/Q \2 r - T l{l + D — 1 ) ~ 2 

’ V 1 - f[r}{df$i mi ) 2 - ^ - VI - f[r] - - p 


(Al) 


where f = r/r h , £ = £/r h , t = r/r h , 4> /m = r h <F irn are dimensionless. 

Performing the following infinitesimal transformation |3S] in the above resultant action: 


f f' = t + e, £->£'= £, 

®l mi (f,£) -A &l mi ( f ',£') = ®lm.i(T,£), 

r{f',£') = r(f + e,£) 


(A2a) 

(A2b) 

(A2c) 


The action in Eq. (|Al| becomes, 


S ~ - y£ / drd£ (r+ ehi + e 2 h 2 /2) 1 


, f h d f e 2 


L df 2 <9 2 / 

h2 dh +hl d^ 


1/2 


(<9f$zmJ 2 




, df L2 d 2 f 

h2 d? +hl d^ 

h 2 % + hf^ 

or or z 


- 1/2 


(d^ lm J 2 - Ki + D~i) _ 

(r + ehi + e 2 /i 2 /2) 


1/2 




(A3) 


where h, = % and h 2 = % 
or or z 
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Using the relation between the Lemaitre coordinates [57] 

£_t= f dr 

J v 71 - /KuO] 


(A4) 


gives the following expression, 


hi 


-x/l 3 /) ^2 


-19/ d|. _ 1 

2 df’ dr' con - T y/r=j 


The Hamiltonian (H) corresponding to the above Lagrangian is 



l .m a 





g? / 2 vr 2 ) 


1(1 +D- 1) 2 
' o X; T 

^1 


(A5) 


(A6) 


where 


/ df e 2 

Si = r + ehi + e 2 h 2 /2, & = W1 - / - efti ^4 - — 


' df ^d 2 f 
h 2 -+~ + ft? 

<9r 1 dr 2 




(A7) 


and n^ mi is the canonical conjugate momenta corresponding to the field x; mi . 

Upon quantization, H; mi and xj mi satisfy the usual canonical commutation relation: 


Xlmi (i,?) jftj'm' = iSu'S mir n>.5 (f - 


(A8) 


Using relations (A5) and expanding the Hamiltonian up to second order in e, we get, 

r 


h - It I. 


~D 


Hr. 


df 


&lr 


(1 + eld 3 - e 2 -ff 4 ) 1/2 f>/» (1 - e-ffj - e 2 ff 2 ) D/2 (1 + eH 3 - e 2 ld 4 ) 1/4 
i(Z + D- 1) 


f 2 (1 — eH 1 -e 2 H 2 ) 


_2 

2 


(A9) 


The Hamiltonian in Eq. (A9|is of the form 


~ i7 0 + el/ + e 2 U 2 


(A10) 


where H 0 is the unperturbed scalar field Hamiltonian in the flat space-time, V, and V 2 are the perturbed parts of the 
Hamiltonian given by; 


//,. 


- £ 


df 


7T 2 + f° 

~lm,i ' ' 


e"/r. 


rD/2 


2 , l(l + D- 1) 2 

H-UA- 


Irrii 


(All) 


U = 




df 


(~Dcjim . + 2 fa' lm .) ( DH 3 aim. - 2 fH 3 a[ m . + 2HfH/ <r, m . - fH' 3 a imi ) 


4r 2 


2i(Z + 14 - I)!?, 2 

H-^- cr, 


^ E 


/»oo 


Irrii 


(A12) 


df 


+ ff 4 ) a& 4 + ( -^ - ™ + DHiH[ - DH 3 H[ + DH' 2 + H 3 H' 3 + \h' a ) a' lmi a lmi 


f 1(1 +D- 1) (3 H\ + 2 H 2 ) tfHl DHh _ D 2 H 3 H( D 2 H 3 H( 1 2 2 _ D 2 H’ 2 

y f 2 4r 2 4f 2 2f 2f 4 1 2f 

-2f " l DH ^ + 16^ ~ IT J 


(A13) 
















































where 


H, = 


vT=7 


H =L d _i 

2 4 f df ’ 


and the redefined field operators are 


1 I / f; — 


= 


Tpr; 


1 9/ 

\/W^ 


iT = 


-1 


4(1 - /) 


a/ 

dr 


1 

2 dr 2 


and Ximi = 


^ Zra» 


( i-/) 1/4 Alm ‘ ' (i - /) 1/4 

such that they satisfy the following canonical commutation relation 

[7T* mi (f, f), CT;/ m / (f', f)] = iS w 6 m . m >.S(f - f') 


(A14) 


(A15) 


(A16) 


The Hamiltonian if in Eq. (AlOl is mapped to a system of N coupled time independent harmonic oscillators (HO) 
with non-periodic boundary conditions. The interaction matrix elements of the Hamiltonian can be found in the 
Ref.[40]. The total internal energy (E) and the entanglement entropy ( S a ) for the ground state of the HO’s is 
computed numerically as a function of e by using central difference scheme. 


Appendix B: Central Difference discretization 


This is one of the effective method for finding the approximate value for derivative of a function in the neighbourhood 
of any discrete point, Xi = Xq + i h,with unit steps of h. The Taylor expansion of the function about the point Xq in 
the forward and backward difference scheme is given respectively by, 


which implies, 


f{x + h) = /( x) + 
f(x — h) = f{x) - M-M. 


, h 2 f"(x) 

2 ! 

, h*f"(x) 
2 ! 


f'(x) 

/"(*) 


/'"(*) 


f{x + h)- f(x - h) + 0 ^ 2 

f(x + h)~ 2 f(x) + f(x - h) + 2 

f(x + 2 h) - 2 f(x + h) + 2f(x - h) + f(x - 2 h) t 

2 h 3 + 1 ’ 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


Appendix C: 


Plots of internal energy and EE as a function of e for different black hole space-times 


In this section of Appendix, we give plots of EE for different black hole space-times; 




Figure 3: Plots of the EE as function of e for the 4-d Schwarzschild black hole with N = 300, n = 50, 80,100, and 130, 
respectively. The blue dots are the numerical data and the red line is the best linear fit to the data. 
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Figure 4: Plots of the EE as function of e for the 4-d Schwarzschild black hole with N = 400, n = 110,140,180, and 200 respectively. 
The blue dots are the numerical data and the red line is the best linear fit to the data. 



Figure 5: Plots of the EE of 4-d R-N black hole in terms of e for different q’s with N = 300, and n = 50. The blue dots are the numerical 
data and the red line is the best linear fit to the data. 
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